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ABSTRACT 

Random matrix models based on an integral over supermatrices are proposed 
as a natural extension of bosonic matrix models. The subtle nature of superspace 
integration allows these models to have very different properties from the analogous 
bosonic models. Two choices of integration slice are investigated. One leads to a 
perturbative structure which is reminiscent of, and perhaps identical to, the usual 
Hermitian matrix models. Another leads to an eigenvalue reduction which can be 
described by a two component plasma in one dimension. A stationary point of the 
model is described. 
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1. Introduction 



Integrals over random matrices have had a variety of physical and mathematical 
applications. They were originally proposed as a statistical model for studying the 
distribution of energy levels of highly excited states of nuclei [1-3]. The matrices 
arc taken to be from some ensemble of diagonalizable matrices, usually Hcrmitian 
or unitary. The pcrturbative expansion of the free energy was found to generate 
Feynman rules which are useful in solving numerous graph-counting problems [4-6]. 
More recently, these properties have been found to make certain scaling limits of 
matrix models useful for studying two-dimensional gravity, or string theory in "less 
than one dimension [7- 11]." 

It is natural to consider the extension of these models to supermatrices, which 
represent the linear transformations of a vector space having both even (bosonic) 
and odd (fermionic) coordinates. All of the algebraic and analytic methods needed 
to define ordinary matrix models have extensions to supermatrices.^ Moreover, 
supermatrices come in two types: c-type, which preserve the coordinate type of a 
vector, and a-type, which interchange fermions and bosons. This suggests that a 
model which includes both could be useful in defining subcritical superstrings. 

This paper will concentrate on c-type matrix models, since these are the most 
closely analogous to the ordinary models. In particular, "almost all" c-type matri- 
ces arc diagonalizable (a-type ones need not be), and diagonalization is a powerful 
tool in rendering bosonic matrix models tractable. Unfortunately, the supcrunitary 
angular integral generally does not decouple, so some nonperturbative methods, 
such as the orthogonal polynomial method, will lose their power. 

Since a c-type supermatrix can be decomposed into a two by two block matrix, 
with bosonic matrices on the diagonal and fermionic ones off it, one may suspect 
that they already contain enough ingredients to define a discretized superstring 
model. However, this seems unlikely to be a proper interpretation, since the po- 
tential contains only bosonic coupling constants, and since the eigenvalues are all 

t A general reference, and the source of all terminology and conventions used here, is ref. [12]. 
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bosonic. If one thinks of c < 1 strings as integrable systems[13], what is wanted 
is a superintegrable system[14], which requires a model with fermionic couphng 
constants to generate both even and odd KdV-type flows. 

One of the greatest distinctions between supermatrix models and ordinary ones 
is that superspace integration is inherently ambiguous. There is no geometrically 
natural reason to restrict the bosonic components of the supermatrix to be pure 

real or complex numbers. In general, they arc simply even elements (c-numbers) in 
a Grassmann algebra, consisting of a body which is an ordinary number, and a soul 
made of even products of anticommuting numbers[12]. To define the integration 
domain requires specifying a slice through soul-space, and different choices really 
should be thought of as different models, with different properties. 

Section 2 of this paper describes the general properties of supermatrix models. 
Section 3 describes the case when the even entries in the matrix are ordinary 
numbers. This choice gives the simplest perturbative expansion, which turns out to 
be very similar (perhaps identical) to that of an ordinary bosonic model. Section 4 
describes integrals over matrices which are constrained to have pure real eigenvalues 
with no soul. (These are called physical supermatrices[12].) When gauge-fixed, 
such a model can be described by a physical system analogous to the Dyson gas of 
bosonic models. In this case, the physical analog is found to be a two-component 
plasma in one dimension. Although the technology for solving such a matrix model 
is presently limited, a saddle point evaluation based on interacting dipoles is used 
to illustrate some features of the model. 

Three very recent supermatrix model references may be of interest. A paper by 
C. Vaz[15] analyzes a special case of the antisymmetric supermatrix model. One 
by L. Alvarez-Gaume and J. Manes, received after the completion of the present 
work, compares the Hermitian matrix and supermatrix models. Finally, G. Gilbert 
and M.J. Perry[16] propose a quenched c-type supermatrix model which appears 
to incorporate genuine supcrsymmctry. 
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2. Supermatrix Integration 



Supermatrix models can be defined by taking any bosonic model and replac- 
ing all of the ingredients by their superspace analogs. In particular, consider a 
Hermitian c-type supermatrix[12] (the only type to be used in this paper) 

where A and C are bosonic and Hermitian, and B is fermionic. The notation 
indicates that i,j,--- denote a general index, while /i,!/, . . . are bosonic and a, (3, . . . 
are fermionic. If (—1)* is +1 for bosonic indices and —1 for fermionic ones, then 
the supertrace is defined to be 

strM = ^(-)*Mii = trA-trC . (2.2) 

i 

For an arbitrary potential V{M), the random Hermitian supermatrix model is 
defined by the partition function 

ZmniP) ^ JdM exp(-/3 stry(M)) , (2.3) 
s 

where M is an (m|n) supermatrix acting on vectors with m commuting and n 
anticommuting components, and 

dM = Yl dA^^ dCaf3 d'B^p . (2.4) 

The measures for A and C are the usual linear Hermitian measures. The integration 
domain S is chosen so that the bodies of A and C range over the ordinary Hermitian 
matices, but the souls, which can contain even products of B components, will 
be specified later. Since the supertrace is indefinite, the integral (2.3) will not 
actually exist for polynomial potentials. However, this is a technical problem that 
can be avoided by considering superunitary matrices instead, or by multiplying the 
potential by i and inserting convergence factors. 
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When S is chosen such that the Gaussian integral (quadratic potential) can 
be evaluated, more general potentials can be handled by a perturbative expansion. 
This is true in particular when A and C are taken to be ordinary soul-less Hermi- 
tian matrices, the case considered in section 3. The present section will consider 
methods which may be applied to more general cases. 

In the bosonic models, the most powerful methods rely on diagonalization to 
obtain an integral over eigenvalues alone, factoring out the degeneracy due to the 
unitary angular integration [3, 5]. Therefore, it is of interest to see what may be 
gained by diagonalizing the supermatrix. A c-type hermitian supermatrix can be 
diagonalized by a superunitary transformation, except in certain singular cases. 
The singular cases occur when one of the eigenvalues of A has the same body as 
an eigenvalue of C. In that limit, a pair of eigenvectors become bodyless, and 
cannot be part of an orthonormal basis, so the diagonalization cannot be carried 
out. However, this occurs only at isolated points in the integral. 

Assume now that 



where U is U{m\n), with a possibly unusual choice of souls. The body of V is 
just the set of (m\n) real diagonal matrices, and all off-diagonal elements vanish. 
Following the classic matrix-model method[3,5], the integrand of (2.3) may be 
multiplied and divided by 




(2.5) 




(2.6) 




U ^ MU and using the invariance of the measure 



dU / rfA A(A)exp(-/?stry(A)) 



(2.7) 



where A is the diagonal matrix of eigenvalues Aj of M. 
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An important difference from the bosonic case is immediately clear. The inte- 
gral over the supergroup does not factorize. The choice of souls for the eigenvalues 
depends on the group element U. If this were a one-dimensional integral, the choice 
of souls would not have mattered, up to a surface term* But for multiple integrals, 
the dependence is nontrivial. Also, the evaluation of (2.6) depends on U, which is 
algebraically U{m\n), but may have unusual soul geometry. 

A case where (2.7) can be evaluated may be obtained by restricting the eigenval- 
ues to be pure real numbers, and U to be the standard [/{mln) with its Haar mea- 
sure. Then the set S consists of what de Witt [12] calls the "physical supermatrices'.' 
This case will be the subject of section 4. It has the drawback that the constraint 
of having real eigenvalues cannot easily be included in (2.3), so that even the Gaus- 
sian integral becomes complicated, and perturbative Feynman rules are difficult to 
develop. 

A possible compromise would be to diagonalize only A and C, by performing 
a superunitary transformation U xV e U (m) xU{n). Then B — > C/Syt. The Ja- 
cobians for A and C are the usual Vandermonde determinants, while the Jacobian 
for B is just 1, so 

fl[{ca-cpfexp{-(3stTV{M)) (2.8) 

where M is now defined with A — diag(a), C — diag(c), and the volume of U (m) x 
U{n) has been divided out. The remaining integral over B can be replaced by an 
integral over a U{m\n) /U{m) x U{n) coset, if desired. It is not clear whether (2.8) 
is a useful supplement to (2.3). This question will not be persued here. 

A simple example may be helpful to clarify the issues in this section. Consider 
the case m — n — 1. Then A and C become real even supernumbers a and c, while 



★ See p. 7 of [12] for a proof. 
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B becomes a complex odd Grassmann number /3. The eigenvalues of M are 

X^a + (a-c)-'^pp*, iJ = c+{a-c)-'^Pl3*, (2.9) 

and is diagonal when the superunitary matrix is 

(1 — haa* a 
-a* 1 + ^aa* 

with a = {a — c)~^(3. 

The partition function will now have the form 



Z= J dadcdpdp*exp{-strV{M)) 

dada* J d\dn{\ - /x)"^ exp(-1/(A) + V{ij)). 



(2.11) 



This is a special case of (2.7) with dU = dada* and A(A,/x) = (A - /x)-^. The 
unimportant diagonal subgroup of t^(l|l) is omitted. Since no even parameters 
are needed in U, there is no ambiguity in expressing the group manifold U in this 
case. The supernumbers a and c can be decomposed into body and soul, which 
have the form 

a + as ^ aB + /(as, cb)PP* 

(2.12) 

c =CB + c's = CB + g{aB, cb)PP*, 

where aB and cb range over the real numbers, and / and g are ordinary real- valued 
functions. This form is unique if the souls depend only on /? and /?*. 

The integral can be reduced to an integral over the bodies by including a 
Jacobian factor 

oo 

Jdadc= j daBdcB[l + {1,1+9,2)1313*]. (2.13) 

— oo 

The presence of such factors can make a direct evaluation of (2.3) complicated for 
large matrices, even for the simplest potentials. It is useful to choose S so that 
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cither (2.3) or (2.7) is as simple as possible. These two cases are the subject of the 
following sections. 

3. Ordinary Hermitian Supermatrix Model 

The partition function (2.3) can be expressed most simply when A and C are 
chosen to be ordinary soul-less Hermitian matrices. This case will be referred 
to as the "ordinary Hermitian supermatrix model'.' No troublesome Jacobian fac- 
tors of the form (2.13) complicate the evaluation, so the Gaussian integral can be 
evaluated exactly. More complicated potentials can be handled perturbatively, by 
developing Feynman rules [5]. The identification of dense Feynman diagrams with 
Riemann surfaces is the reason for the relevance of a scaling limit of ordinary ma- 
trix models to two-dimensional gravity [7-9]. Therefore, it should be expected that 
the connection to gravity can be made most directly for the ordinary Hermitian 
supermatrix model. 

Consider a quadratic potential |strM^. The Feynman rules and Wick's the- 
orem can be derived as in the usual case, by calculating the expectation value of 
a source exp(str ( JM)), and varying with respect to J to obtain the expectation 
values needed for a perturbative expansion. The propagator is simply 

(MyM,*,) = i-ySikSji. (3.1) 

Wick contractions must always be carried out after permuting the matrices so that 
they are adjacent, keeping track of the signs (— )(*+j)('=+0 introduced by commuting 
My past M]fi. Symmetry factors can be calculated as in the usual matrix models, 
following the rules of ref . [5] . 

V(M) = istr M^ + g str , (3.2) 
the lowest-order contribution to the free energy comes from the two-loop graphs 
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Figure 1. The three lowest-order graphs for the quartic model, 
in figure 1, which represent (strM^). Evaluating them gives 

g (str =g J2 i-y [{MiMj) (MklMli) + (M.^-Mf,) {Mj^M^i) 

i,j,k,l 
i,j,k,l 

—g [(m — n)^ + (m — n) + (m — n)^] . 

(3.3) 

The three contributions to (3.3) come respectively from the three graphs in figure 
1. Graphs (a) and (c) are planar, while (b) is nonplanar. It was possible to 
rescale g so that (3.3) is identical to the standard Hermitian matrix result[5]. If 
g — g{m — n)~^, then (3.3) may be written 

g{m - n)~^ (strM^) = g[2{m - nf + l] (3.4) 

The power of {m — n)^ is (1 — G), where G is the genus of the surface identified 
with the graph by filling in faces in the usual manner. 

This provides evidence that the supermatrix model will have a 2d gravity in- 
terpretation, when m,n — > oo with n/m fixed. A computation of all three-loop 
connected graphs containing two quartic vertices provides further evidence for this. 
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Figure 2. The inequivalent connected three- loop graphs with two quar- 
tic vertices. 

The total contribution of these graphs, shown in figure 2, to the free energy is 

— — ((str M^f) = g\l8(m - nf + 30l , (3.5) 
2!(m-n)2 ' /conn. yi V I J > \ ) 

which is identical to the ordinary Hermitian matrix result [5], with N replaced by 
m — n. In fact, it is clear that in an arbitrarily complex diagram, each closed 
index loop (face in a generalized triangulation of the surface) will give a factor 
of m ± n. If the pattern suggested by low order computations persists, then only 
m — n will occur, and the model will be precisely equivalent to the \m — n\ x \m — n\ 
Hermitian matrix model. (Ref. [17] provides some further evidence that this is the 
case.) However, if any factors of m + n occur, the topological expansion will still 
go through, but the weighting factors for various surfaces will be changed. More 
powerful combinatorial methods are needed to settle this question. 

If the identification with the |m — n\ Hermitian model persists, then this may 
suggest a connection with the bosonic sector of a subcritical superstring theory. 
Ordinary superstrings have a bosonic Nevcu-Schwarz sector which by itself is closely 
analogous to the purely bosonic string, in the sense that is can contain the same 
massless spacetime fields (gravity, etc.) with the same low-energy effective action. 
The fermions play a purely internal role in this sector. This seems analogous to 
the role played by the fermions in a c-type supermatrix. They play an internal 
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role, but in a "spacetime" sense, the model behaves as a bosonic one. Adjoining 
a-type supermatrices could provide an analog of a Ramond sector, containing true 
fermions. Quenching, or constraining some elements of a supermatrix element, 
may also lead to a supersymmetric theory [16]. 

4. Physical Supermatrix Model 

Having considered the case where the even entries in the supermatrix are pure 
numbers, it is now natural to turn to the other relatively simple case, where the 
eigenvalues of the supermatrix are pure real numbers. Such Hermitian superma- 
trices may be thought of as observables in a quantum-mechanical system defined 
on a supermanifold* 

Incorporating the constraint that M have real eigenvalues directly into (2.3) 
is complicated. A Jacobian factor will be needed, which is difficult to compute 
in general. However, working with the eigenvalue reduction (2.7) is simple in this 
case. The eigenvalues Aj are integrated over the ordinary real numbers, and the 
supergroup integral over U{m\n) factorizes. In fact, that integral vanishes, as can 
be seen from (2.11) for the (1|1) supermatrix example, or in general from the fact 
that (2.3) vanishes when dXi/dB^a — 0. In any case, the supergroup integral 
will be dealt with by fixing the gauge, and simply dropping it. This gives the 
closest analog of the ordinary matrix model eigenvalue reduction. The remaining 
eigenvalue partition function will be essentially a vector model whose form and 
symmetries are inherited from the underlying supermatrix model. 

The only ingredient still needed in (2.7) is the determinant A. This can be 
calculated by representing U in (2.6) as e-^ with A anti-Hermitian. Then, if A is 
the diagonalized matrix, 

(C/AC/t).. = Ay(A,--A,) + ... (4.1) 

and, keeping track of the Grassmann character of the delta functions in (2.6), the 
★ Such systems are described in chapter 5.3 of ref. [12]. 
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determinant is found to be 

A(A) = n (-^M - ^-)' n - ^/^)' IK^M - (4.2) 

^l<^' a<l3 IJ^a 

in the index conventions of (2.1). The poles occur at supermatrices which tech- 
nically are not diagonalizable for the reasons noted in section 2. The gauge- fixed 
physical supematrix model is then defined by the equation 

/ rn n 
H dXiA{X) exp -/? (J2 ^(Am) - E ^(^«)) ■ (4-3) 
i /i=l a=l 

This model is clearly inequivalent to the usual hermitian model, due to the 
presence of a denominator in (4.2). The ordinary matrix models were introduced 
as a model for describing the eigenvalue distribution of a random physical operator. 
In that case, Wigner[l] observed that the eigenvalues repel each other, and Dyson[2] 
noted that they can be interpreted as a gas of charged particles with = in 
one dimension, at temperature In the limit N, f3 ^ oo with f3 < N, the 

gas freezes into a crystal with a charge distribution determined by the external 
potential V. The critical case takes P/N — > 1, so that fluctuations about the 
crystal become important, and in the perturbative expansion, higher genus graphs 
contribute comparably to planar ones [7-9]. 

The physical analog for the supermatrix model is a two-component plasma in 
one dimension, consisting of equal and opposite charges ±/3^^/^ at temperature 
in an external potential V. The Ay^^ may be considered to be m positive 
charges, while the Xa are n negative charges. At very low temperatures, the positive 
charges will pair with a negative partner, leading to an effective theory of n dipoles 
interacting with m — n positive charges when m > n. 

Assume m > n and introduce variables 

Xi — ^(Xi + Xffi+i), Ui — ^(Xi — Xfn+i) , Ha — Xn+a (4-4) 

ior i — 1, . . . ,n, and a — 1, . . . ,m — n. In the limit where the first n positive and 
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negative charges pair, Xi is a center of mass coordinate for the ith. dipole, Ui is its 
charge displacement, and ya labels the m — n left-over positive charges. In terms 
of these variables, (4.3) may be rewritten 



7 _ on 



» n , , m—n 



n 

i,a 



1 « o=l 

l+Ui{Xi - ya) 



1- {Ui- UjY{Xi - Xj] 



-ln2 



1- {Ui + UjY{Xi - Xj) 2 
a,b 

X exp -/? I ^(^«) + Yl + ~ ~ ^*)] 

La i 



1 - Wi(xj - ya) 



(4.5) 



Clearly, Z^n.^^ is dominated by configurations where -Uj — ^ 0, and by configura- 
tions equivalent to this up to index permutations. These configurations actually 
cause Z to diverge, which is a price for simply dropping the U{m\n) integral. 
However, once the internal energy of the dipoles is regulated by cutting off the 
displacements, no further divergences occur, provided the potential is well-chosen. 

For simplicity, further attention will be restricted to the neutral plasma, with 
m = n. In the dipole approximation, cutting off the displacements at 6i gives 



oo 



-oo Si 



(4.6) 



X 



Y exp -/5 <{ ej [V{xi + Ui) - V{xi - Ui)] + ^ ^i^j^ij 

i i^j 



ei=±l 



where Ui are now considered to be small and positive and the dipole interaction is 



^ —1 
Vij == — — tanh 

P 



2uiUi 



— — — 



(Xi X 



3) 



UiUj 



(5 (xi •^j)' 



(4.7) 



To lowest order in the interaction is attractive, but the exact form shows that 
higher order repulsive effects will prevent a singularity in Vij for The 
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dipole interaction with the external potential is eiUiV'{xi), where V' is the deriva- 
tive of V with respect to x. Therefore, positively polarized dipoles are drawn to 
the minima of V' while negative ones are drawn to its maxima. 

It is informative to choose a simple potential, and find the behavior of the free 
energy near a saddle point. Consider the example of a cubic potential 



In this case, the true minima of the potential should come from unstable configura- 
tions where negative dipoles lump together and run off to infinity, or positive ones 



tamed by modifying the potential at large (a situation familiar in matrix mod- 
els). The latter case is important, but requires the inclusion of the higher-order 
repulsive effects in Vij for meaningful results. 

Therefore, an unstable stationary point, where negative dipoles are repelled 
from the minimum of V' but attracted to each other, will be considered for illus- 
trative purposes. (This is opposite the usual situation for matrix models.) It is a 
simple enough case to allow a direct comparison with the saddle point analysis of 
the standard matrix model[4] , and it is useful to know exactly how the supermatrix 
analysis differs. 

In the n, /? ^ oo limit, it is useful to define a continuum variable x{t) — n~^l'^Xi 
with t — i/n. Near the saddle point, the dipole approximation should be valid, 
and we will work with the effective dipole energy 



l3V{x) = \x^ + \(5-^l'^gx'^. 



(4.8) 



all condense at x 



^. The former case will be neglected, since it can be 
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and cutoff in (4.6), and for the cubic potential, 

/^\ 1/2 

(3V'{x{t)) = x{t) + gx\t), g = ^-j g. (4.10) 

Introducing a dipole density P{x) — p{t)dt/dx and dropping the dipole self- 
energy term gives 

n^pE^2 J dxP{x){x + gx^) -S j j dxdx' (4.11) 

The saddle point equation SE/Sp{t) = is satisfied when dx{SE/SP{x)) — 0. For 
(4.11) this imphes 



X 



1 16 /■ , , Fix') ,^ , 

+ 7y^ = -^fdx' ^ 4.12 

^9 9 J [x- x'y 



using the same principle part prescription as ref. [4] to handle the coincident point. 
(Recall that (4.6) shows that the coincident limit is not truly singular.) 

Since the saddle point which will be found is the unstable one with negative 
dipoles clustered about the minimum of V' (using some foresight), it is convenient 
to define the integrated dipole density as 

J dxP{x) = -po- (4.13) 

The quantity po is determined by how n~^/'^6i is tuned as n — > cx), 5j — > 0, and 
should be thought of as a new parameter of the model. It is also convenient to 
introduce a new variable 

w ^ k(x + , (4.14) 



2? 

centered at the minimum of V', with scale k,~^ — 32pq/^ chosen for later simplicity. 
Then the normalized negative dipole density p defined by 

P(x) = -^'KK,pop{w) (4-15) 

has total integral +2/7r. (This integral was fixed for later convenience.) 
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In the new variables, the stationary condition (4.12) becomes 



w 



= T / d^'r^^^ > kl < «■ (4-16) 



Due to the symmetry of the problem, p{w) will be an even function, with support in 
an interval [—a, a] to be determined. The density does not have to be positive, since 
dipoles can have two polarizations, but the integral must be 2/7r by construction. 
A solution can be found following ref. [4]'s treatment of the ordinary quartic matrix 
model, by introducing an analytic function 

F(v,) = lLv,'jJ^^, (4.17) 
4 J {w-w'Y 

—a 

which is analytic for complex w with cut [—a, a], behaves as \w~^ for \w\ — > oo, 

and satisfies 

F{w ±ie)^w^ Imp'^w) (4.18) 

for w in [—a, a]. 
The solution is 

F{w) = ty - yJiiP- — w~'^, (4.19) 

which requires a — 1 and 

P 



"{w) = — Vl-«;4 (4.20) 



\w\ 



for \w\ < 1. Integrating once gives 

p'{w) = (^Vl - w4 - tanh"^ \/l - sgn w. (4.21) 
The logarithmic singularity at w = is integrable, and the resulting density will 
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Figure 3. The normalized negative dipole density p{w) as a func- 
tion of the coordinate w centered at the minimum of the 
quadratic potential V'. 

be finite. The integration constants are fixed by adding a term 

-h{l - \w\) (4.22) 

to p{w), with coefficient h chosen such that the integrated density is 2/7r, as re- 
quired by the defining conditions. Eqn. (4.21) may be integrated numerically, 
leading to /i = .147, for which the final solution p{w) is shown in figure 3. 

The energy for this configuration is found by substituting p into (4.11). Using 
some integration by parts, (4.20), (4.21), and (4.22) suffice to obtain ^E. The 
double integral can be done with the help of (4.12) as in ref. [4], leading to 

n2/5E = |^-% + i7r^i/V'{^(«) - 3 + /^log(l - «)} (4-23) 
where = and 

I{a) = - / -VT^ log(l -az)^^y , . (4.24) 

{ ^ 2 ^r(^)n(n+l) ^ ^ 

The surface interpretation of a matrix model depends only on the non-analytic 
behavior of the partition function. The expression (4.23) has non-analytic behavior 
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when 



) 



3/2 



1. 



(4.25) 



This can occur for n//? ^ 1 if g^pQ — 2 ^. 

Since this fixed point is unstable, it is not clear that it has any application 
in 2c? gravity. (However, inverted potentials sometimes come up in scaling limits 
of matrix models [9], so instabilities are not always as bad as they seem.) This is 
the only kind of stationary point that can be found within the limits of the dipole 
approximation of the one-dimensional plasma. However, the interaction potential 
Vij contains higher-order repulsive effects, so stable fixed points may eventually 
be found. A more exact treatment would be desirable, but this is more difficult 
than for the usual matrix models, where the orthogonal polynomial method yields 
relatively easy exact results. 



The motivation for introducing supermatrix models is their apparent mathe- 
matical similarity to ordinary matrix models. This suggests that they could provide 
a way to incorporate fermions in a matrix model without sacrificing integrability. 

However, when these models are investigated in detail, it becomes clear that 
the nature of superspace gives them very different properties from ordinary matrix 
models, and that certain useful tools, especially the orthogonal polynomial method, 
are not easily applied. The source of the most unique features of supermatrix 
models is the inherent ambiguity in integrating over even Grassmann variables. 

The most obvious definition of the supermatrix integral gives a perturbative 
expansion which implies that a scaling limit of the model will have a 2d gravity 
interpretation (but not supergravity, since the models considered have only bosonic 
coupling constants). The fact that the superunitary integral does not decouple in 
this case means that the eigenvalue reduction an orthogonal polynomial method 
are not readily applicable, although eventually suitable extensions may be found. 



5. Conclusions 
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The leading terms in the perturbative expansion of the quartic (m|n) Hermitian 
supermatrix model are identical to those of the (m — n) bosonic Hermitian matrix 
model. This may suggest a connection to the "spacetime" bosonic sector of a 
subcritical superstring, to which fermionic partners must be adjoined to complete 
the theory. The fact that supermatrices contain odd Grassmann variables is not 
in itself reason to believe that they would constitute a supersymmetric theory by 
themself. Perhaps the appropriate partners are supermatrices with the bosonic 
and fermionic entries enterchanged. 

Since the supermatrix integrand depends only on invariants, it will actually 
vanish whenever the superunitary integral does decouple. Nevertheless, it is useful 
to define such models by gauge- fixing, since they provide the closest super-analog of 
Wigner's model for the eigenvalue distribution of a physical operator. The result is 
that the eigenvalues of a physical super-operator separate into two classes. Within 
each class, the eigenvalues repel, as in Wigner's model. However, eigenvalues in 
each class attract those of the other. If any real system has such properties, perhaps 
supermatrices could provide a useful model for some of its statistical properties, as 
matrix models did with some degree of success[18] for the highly-excited states of 
nuclei. 

In analogy with the Dyson gas interpretation of matrix models, the eigenvalues 
of a supermatrix may be though of in terms of a two-component plasma in one 
dimension. This physical analog could provide useful intuition for solving the 
matrix model. Conversely, the matrix model, when sufficiently developed, could 
provide a new source of information about such plasmas. Electron-hole plasmas 
in optically-excited semiconductors are a well-studied system[19] which may be 
somewhat analogous. 

The analysis presented here should be considered preliminary. Assuming that 
opposite charges condense into pairs at low temperatures, an effective dipole gas 
was constructed, and one of the stationary points for a cubic potential was ana- 
lyzed. The unstable fixed point analyzed was chosen because it is the simplest. If 
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a way is found to deal with the higher-order repushve effects present in (4.7), it 
may be possible to analyze a stable fixed point instead. It may be necessary to 
extend the analysis to include effective neutral particles containing more than two 
charges, i.e. higher multipoles. Developing a more exact treatment of the plasma 
could reveal some interesting physics, and perhaps new 2d gravity models. 
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Figure Captions 

Figure 1. The three lowest-order graphs for the quartic model. 

Figure 2. The inequivalent connected three- loop graphs with two quartic vertices. 

Figure 3. The normalized negative dipole density p{w) as a function of the 
coordinate w centered at the minimum of the quadratic potential V'. 
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